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g : 1 Introduction 



The investigation of non-perturbative properties of supersymmetric (SUSY) gauge 
theories is presently a rather theoretical subject: it is not known whether (broken) 
supersymmetry is realized in nature or not. Moreover, the simple ("minimal") su- 
persymmetric extensions of the standard model do not involve strong interactions 
near or above the scale of supersymmetry breaking, hence the knowledge of the non- 
perturbative supersymmetric dynamics is not directly required. In spite of this there 
is a continuing interest in studying strongly interacting quantum field theories with 
supersymmetry. (For an early review of the subject see ref. 

The motivation to investigate non-perturbative features of supersymmetric gauge 
theories is partly coming from the desire to understand relativistic quantum field 
theories better in general: the supersymmetric points in the parameter space of all 
quantum field theories are very special since they correspond to situations of a high 
degree of symmetry. As recent work of Seiberg and Witten and other related 
papers showed, there is a possibility to approach non-perturbative questions in four 
dimensional quantum field theories by starting from exact solutions in some highly 
symmetric points and treat the symmetry breaking as a small perturbation. Beyond 
this, the knowledge of non-perturbative dynamics in supersymmetric quantum field 
theories can also be helpful in understanding the greatest puzzle of the standard 
model, with or without supersymmetric extensions, namely the existence of a large 
number of seemingly free parameters. As we know from QCD, strong interactions 
in non-abelian gauge theories are capable to reproduce from a small number of in- 
put parameters a large number of dynamically generated parameters for quantities 
characterizing bound states. This is a possible solution also for the parameters of 
the standard model if new strong interactions are active beyond the electroweak 
symmetry breaking scale. 
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1.1 N = 1 SYM 



The simplest supersymmetric gauge theory is the supersymmetric extension of Yang- 
Mills theory. The action of Yang-Mills theory with N = 1 supersymmetry is conven- 
tionally given as 

' d 4 xd 2 6Tr(W a W a ) 



= J d 4 xTr {-£ V" + % -F^ - iXa^D^X) + i{D^X + -D 2 j , (1) 

where the first line is written in terms of the spinorial field strength superfield 
W(x, 8, 9) a which depends on the four-coordinate x and the anticommuting Weyl- 
spinor variables 6 a ,6a (ck, d = 1,2). After performing the Grassmannian integration 
on 9, one obtains the second form in terms of the component fields. The field strength 
tensor F pv and its dual are defined, as usual, as 

F^x) = -igF; u {x)T r , F, u ee \t, upa F^ . (2) 

A, A represent a Majorana fermion field in the adjoint representation and D is an 
auxiliary field. 

The action in (|l|) includes a B-term, therefore it is natural to introduce the 
complex coupling 

6 Ani 
2^ + J 2 

and then, with arbitrary 0, the N = 1 SYM action becomes: 



C Hill I . . 

T = - + ( 3 ) 



^{r J d 4 x d 2 0Tr(W a W a )} 



— I d 4 x r l i" 

g 2 



-F^F^ - iXa^D^X) + ((D^X + D 2 



+ 7^7 / (I'.rTv 



16tt 2 y ~ ~ - r J (4) 

The remarkable feature of the action in (H) is that, after performing the trivial 
Gaussian integration over the auxiliary field D, it is nothing else than an ordinary 
Yang-Mills action with a massless Majorana fermion in the adjoint representation. 
This shows that this theory is "automatically" supersymmetric. Introducing a non- 
zero gaugino mass rrig breaks supersymmetry "softly" . Such a mass term is: 

mg(X a X a + X & X & ) = m s (M) . (5) 

Here in the first form the Majorana- Weyl components A, A are used, in the second 
form the Dirac-Majorana field 

The Yang-Mills theory of a Majorana fermion in the adjoint representation is, in 
a general sense, similar to QCD: besides the special Majorana- feature the only differ- 
ence is that the fermion is in the adjoint representation and not in the fundamental 
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one. As in QCD, a central feature of low-energy dynamics is the realization of the 
global chiral symmetry. As there is only a single Majorana adjoint "flavour", the 
global chiral symmetry of N — 1 SYM is U(l)\, which coincides with the so called 
R- symmetry generated by the transformations 

0« = e*0 a , O'a = z~ lip h ■ (6) 

This is equivalent to 

X' a = e iip X a , A' d = e~ ilp K , = e^ 6 * . (7) 

The [7(1 ^-symmetry is anomalous: for the corresponding axial current = 
^7^75^5 i n case °f SU(N C ) gauge group with coupling g, we have 



&J = ^e^FLFL . (8) 

327T 



However, the anomaly leaves a Z 2 jv c subgroup of ?7(1)a unbroken. This can be seen, 
for instance, by noting that the transformations 

$ e -^75^ ^ vj> (Q) 

are equivalent to 



nig — >■ m^e 



" 2 ^ 75 , 6 5 yA/ -> Qsym - 2N cV > , (10) 



where Qsym is the 0-parameter of gauge dynamics. Since Qsym is periodic with 
period 2n, for m- g = the U(l)\ symmetry is unbroken if 

<p = ( p k = —, (fc = 0,l,...,2iV c -l) . (11) 

For this statement it is essential that the topological charge is integer. 

The discrete global chiral symmetry Z 2 n c is expected to be spontaneously broken 
by the non-zero gaugino condensate (XX) / to Z 2 defined by {(po,<fN c } (note 
that A — > —A is a rotation). The consequence of this spontaneous chiral symmetry 
breaking pattern is the existence of a first order phase transition at zero gaugino 
mass nig = 0. For instance, in case of A^ c = 2 there exist two degenerate ground 
states with opposite signs of the gaugino condensate. The symmetry breaking is 
linear in nig, therefore the two ground states are exchanged at m g = and there is 
a first order phase transition. 

The non-perturbative features of the SYM theory can be investigated in a lattice 
formulation. As always, the lattice action is not unique (see section EO). A possible 



formulation was given by Curci and Veneziano |y based on the well known lattice 
formulation of QCD introduced by Wilson. In the lattice action the bare gauge cou- 
pling (of SU(N C )) is convetionally represented by (3 = 2N c /g 2 and the bare gaugino 
mass by the hopping parameter K. In the plane of (/3, K) there is a critical line cor- 
responding to zero gaugino mass and the expected phase structure is the one shown 
in figure [I]. 
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Figure 1: Expected phase structure of SU(2) Yang-Mills theory with 
a Majorana fermion in adjoint representation in the (/3, -fT)-plane. The 
dashed-dotted line is a first order phase transition at zero gluino mass, 
where supersymmetry is expected. 



1.2 N = 2 SYM 

The SYM theory with N = 2 extended supersymmetry is a highly constrained theory 
which has, however, more structure than the relatively simple N — 1 case discussed 
above. In particular, besides the N — 1 "vector superfield" containing the gauge 
boson and gaugino (A^X), it also involves an N = 1 "chiral superfield" (0, A') in 
the adjoint representation which consists of the complex scalar = A + %B and the 
Majorana fermion A'. The Majorana pair (A, A') can be combined to a Dirac-fermion 
ip and then the vector-like (non-chiral) nature of this theory can be made explicit. 

The Euclidean action of N = 2 SYM theory in component notation, for simplicity 
in case of an SU(2) gauge group, is the following: 

Sp 2 M = J d 4 x {\f; u {x)F;M + \{D^{x)){D,A r {x)) + \{D ,W \x))(D ,W \x)) 
+tf\x) ltl D^ r (x) + ige rs $ r (x)[A s (x) + i lb B s (x)}^ 1 (x) 

+ ^[A r (x)A r (x)B s (x)B s (x) - A r (x)B r (x)A s (x)B s (x)] } . (12) 

This is a massless adjoint Higgs- Yukawa model with special Yukawa- and quartic 
couplings given in terms of the gauge coupling. 

In iV = 2 SUSY possible couplings are so strongly constrained by the symme- 
try that only gauge couplings are allowed. Another important feature is that the 
symmetry also implies that the matter field content is always vector-like. Therefore 



(AX)<0 



(AA)>0 - - ~-~ b ^les, 
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N = 2 SUSY theories are always non-chiral and hence well suited for lattice sudies. 
A lattice formulation of iV = 2 SYM based on the Wilsonian formulation of QCD 
has been investigated in Q. 

The main new feature of N = 2 SYM compared to N = 1 SYM is that it also 
contains scalar fields, hence there is the possibility of Higgs mechanism. Let us 
here only consider the simplest case of an SU{2) gauge group. In the Higgs phase 
the vacuum expectation value of the scalar field is non-zero. In terms of the real 
components <j>(x) = A(x) + iB(x) we have 

(A r (x)) ? , (B r (x)) ^ • (13) 

This implies the Higgs mechanism breaking of SU(2) — > U(l), similarly to the 
Georgi-Glashow model. Due to the Higgs mechanism the "charged" gauge bosons 
become heavy. The low-energy effective theory is iV = 2 SYM with U(l) gauge 
group. 

Seiber and Witten proved that extended SUSY and asymptotic freedom can be 
exploited to determine exactly the low-energy effective action in the Higgs phase, if 
the vacuum expectation values are large. At strong couplings there are two singular- 
ities of the effective action corresponding to light monopoles and dyons, respectively. 

The expectation value of the complex scalar field (0) parametrizes the moduli 
space of zero-energy degenerate vacua. The degeneracy is a consequence of N = 2 
supersymmetry. This phenomenon is usually referred to as the existence of flat 
directions: the potential identically vanishes for B r = cA r . In the present case the 
moduli space is a non-compact manifold with two parameters. The presence of non- 
compact flat directions requires the breaking of supersymmetry for the definition 
of the path integral over the scalar fields: otherwise the path integral would be 
divergent. It is also plausible that soft breaking with mass terms is not enough 
in the Higgs phase, where the mass-squared terms in the potential are negative. 
Therefore small hard breaking by dimensionless couplings is also required. 

These arguments are quite general. In the special case of N = 2 SYM the general 
renormalizable scalar potential with the given set of scalar fields is 

V(A,B) = \m\A r A r + ^m 2 B B r B r 

+ X A {A r A r f + X B (B r B r ) 2 + X [AB] A r A r B s B s - X (AB) (A r B r ) 2 . (14) 
N=2 supersymmetry is at the point of parameter space where 

m A = mB = \ A = \ B = , 2X [AB] = 2\ {AB) = g 2 . (15) 

In order that the path integral over the scalar fields be convergent, the quartic 
couplings have to fulfil the following conditions: 

X A > AND \ B > AND 
{ X [AB] > max(0, X(ab)) OR AX A X B > max[Aj^ B] , (X [AB] - X {AB) ) 2 ] } . (16) 
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This is in conflict with the super symmetry conditions. 

The consequence of the conflict between supersymmetry and the convergence of 
the path integral over the scalar fields is that in a path integral formulation of the 
quantized theory the supersymmetry has to be broken. On the lattice this means 
that supersymmetry is broken as long as the lattice spacing is non-zero and can only 
be restored in the continuum limit a — > 0. 

The tuning to the supersymmetric point in N = 2 SYM can be studied, for 
instance, in lattice perturbation theory M. It can be shown that the compact flat 
direction is reproduced for a — > on a specific phase transition where three different 
kinds of Higgs-phases meet. The emergence of the non-compact flat direction is a 
result of cancelling of quantum correction contributions from scalars and fermions. 
This is similar to the situation which occurs if the so called vacuum stability boundary 
in Higgs- Yukawa models reaches zero fermion mass. (For a lattice investigation of 
the vacuum stability bound see, for instance, ref. ||.) 

2 Non-perturbative predictions for N = 1 SYM 

In analogy with QCD, one expects that the spectrum of the SYM model consists of 
colourless bound states formed out of the fundamental excitations, namely gluons 
and gluinos. (In this context we shall use the name "gluino" instead of the more 
general term "gaugino".) In the supersymmetric point at zero gluino mass these 
bound states should be organized in supersymmetry multiplets, according to the 
representations of the SUSY extension of the Poincare algebra. For the description 
of lowest energy bound states one can use an effective field theory in terms of suitably 
chosen colourless composite operators. 

For N — 1 SYM the effective action was constructed by Veneziano and Yankielow- 
icz ( VY) . The composite operator appearing in the VY effective action is a chiral 
supermultiplet 5* containing as component fields the expressions for the anomalies 

S = A(y) + V29V(y) + 6 2 F(y) , (17) 
where, for instance, the scalar component is proportional to the gluino bilinear 

A oc X a X a . (18) 

The other components contain gluino-gluino and gluino-gluon combinations. There- 
fore, as far as a constituent picture is applicable to the bound states formed by 
strong interactions, the particle content of the lowest supersymmetry multiplet is: a 
pseudoscalar gluino-gluino bound state, a Majorana spinor gluon-gluino bound state 
and a scalar gluino-gluino bound state. In terms of S the VY effective action has 
the form 

S VY = - [ d 4 xd 2 6d 2 6(S^S) 1/3 + ^[1 d*xd 2 6(S log ^- -S) + h.c] . (19) 
oi J J A 
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Here a and 7 are positive constants and A is the usual mass parameter for the 
asymptotically free coupling defined at scale \i: 

A = ^-i/VoM? , ft = ^ . (20) 

As usual, Pq denotes the first coefficient of the /5-function and we consider here the 
gauge group SU(N C ). 

The effective action in ( (L9|) incorporates the breaking of the discrete Z2N C chiral 
symmetry by the gluino condensate 

{XX} = CA 3 e 2nik/Nc . (21) 

The phase factor depending on the integer k refers to the different ground states 
defined in (|TT|). The proportionality factor C depends, of course, on the renormal- 
ization scheme belonging to A. Instanton calculations and other reasonings imply 
that we have C = 32tt 2 in the dimensional reduction scheme A = A^ [[|. 

The main assumption needed to derive the VY effective action (|20|) is the choice 
of the chiral superfield S as the dominant degree of freedom of low energy dynamics. 
Making a more general ansatz also containing gluon-gluon composites leads to a 
generalization and to two mixed supermultiplets in the low energy spectrum ||. 
Even if S is accepted as the dominant variable, one can argue about the existence of 
a chirally symmetric ground state, in addition to the N c ground states with broken 
chiral symmetry given by the integer k < N c in (|2~T|) []10[. 



An interesting question is how the spectrum of glueballs, gluinoballs and gluino- 
glueballs is influenced by the soft supersymmetry breaking due to a non-zero gluino 
mass rrig 7^ 0. For small m g it is possible to derive the coefficients of the terms 
linear in m g in the mass formulas |Tl| . (Note, however, that the two papers in this 
reference arrive to different results.) 

A general consequence of the chiral symmetry breaking is the existence of a first 
order phase transition at m g = 0. At this point the different ground states in ( pip 
are degenerate and a coexistence of the corresponding phases is possible. In a mixed 
phase situation, as usual at first order phase transitions, the different phases are 
separated by "bubble wall" interfaces. The interface tension of the walls can be 
exactly derived from the central extension of the N — 1 SUSY algebra [|EJ. The 
result is that the energy density of the interface wall is related to the jump of the 
gluino condensate by 

6=^|(AA) 1 -(AA) 2 | . (22) 

Combining this with eq. ( |2"T| ) implies that the dimensionless ratio e/| (AA) | is predicted 
independently of the renormalization scheme. 

In order to compare the predictions (^Tj) and (22) to the results of lattice Monte 



Carlo simulations it is convenient to switch to the lattice A-parameter Alat- First 
one can use || 

WAMs = exp{-l/18} (23) 
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and then for the Curci-Veneziano lattice action |p~3 



A m /A L AT = expl- — 



N n 



iV 

ft = — ^(H - 2n a ) , P = 0.0849780(1) , P 3 = 0.0066960(1) . (24) 

Here n a is the number of Majorana fermions in the adjoint representation, that is 
for SYM we set n a — 1. 

For transforming fl2"T| ) and (E3) to lattice units we need, in fact, the value of 
ah-LAT at the particular values of interest of the lattice bare parameters /3, K. Before 
performing the lattice simulations this is, of course, not known. An order of mag- 
nitude estimate can be obtained from pure gauge theory (at K = 0) by noting that 
both for N c = 2 and N c = 3 we have for the lowest gluball mass M 0] 

aM 

aA LAT ^ — . (25) 

Assuming this approximate relation also at the critical line for zero gluino mass, 
we can use eqs. fl2"ID-([25]) for estimating orders of magnitudes. In the region where 
aM = 0(1) we obtain 

a 3 |(AA)! - (AA) 2 | = 0(1) , a 3 e = 0(10^) . (26) 

As these numbers show, the predicted first order phase transition is, in fact, strong 
enough for a relatively easy observation in lattice simulations. 



3 Numerical Monte Carlo simulations 

The lattice Monte Carlo simulations of quantum field theories are performed in 
Euclidean space-time. For SYM, and more generally for a Yang-Mills theory of 
Majorana fermions in the adjoint representation ("gaugino" or in the context of 
strong interactions "gluino") with arbitrary mass we need first of all the definition 
of Majorana fermions in Euclidean space-time. 

In the literature one may sometimes find the statement that there are no Eu- 
clidean Majorana spinors (see, for instance, |15|). This is only true as long as one 
is concentrating on the hermiticity properties of fields, as in Minkowski space. The 
definition required for an Euclidean path integral can be based on the appropriate 



analytic continuation of expectation values |TJ|. The essential point is that for Ma- 
jorana fermions the Grassmann variables \I/ and \& are not independent, as is the 
case for Dirac fermions, but are related by 



* = * J C , (27) 

with C the charge conjugation Dirac matrix. In fact, starting from an Euclidean 
Dirac fermion field represented by the pair if), if> one can define two Majorana fermion 
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fields satisfying (p7| ) by 



tf« = -L(Vi + C^ T ) , M> (2) s -^(-Vi + C/) . (28) 



Using also the inverse relations 



V = + ^ (2) ) , ^ s Cij? = -±=(*« - zM/^) , (29) 

one can easily relate expectation values of Majorana and Dirac fermion fields 



3.1 Lattice actions and algorithms 

Following Curci and Veneziano ||, we can take for the fermionic part of the SYM 
action the well known Wilson formulation. If the Grassmanian fermion fields in the 
adjoint representation are denoted by ip r x and ip x , with r being the adjoint represen- 
tation index (r = 1, ..,N% — 1 for SU(iV c ) ), then the fermionic part of the lattice 
action is: 

s f = E - k E pw A u w (i + -yM + - 7*)V£+ A ] } • (3°) 

a; /i=l 

Here K is the hopping parameter, the irrelevant Wilson parameter removing the 
fermion doublers in the continuum limit is fixed to r = 1, and the matrix for the 
gauge-field link in the adjoint representation is defined as 

V rs , XfX = V rs ^[U] = 2Tr(UlT r U xll T s ) = V r * s ^ = . (31) 

The generators T r = |A r satisfy the usual normalization Tr (A r A s ) = |. In case of 
SU(2) (N c = 2) we have T r = |r r with the isospin Pauli-matrices r r . The normal- 
ization of the fermion fields in (|3~0D is the usual one for numerical simulations. The 
full lattice action is the sum of the pure gauge part and fermionic part: 

S = S 9 + S f . (32) 

Here the standard Wilson action for the SU(iV c ) gauge field S g is a sum over the 
plaquettes 

S 9 = P52(l-±-ReTrU pl ) , (33) 
pi v iVc J 

with the bare gauge coupling given by (3 = 2N c /g 2 . 

Using the relations in fl29|) one can decompose Sf as a sum over the two Majorana 
components: 

s f = E = \ E E *r^,^ i)u > ( 34 ) 

xu,yv j=\ xu,yv 
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where the fermion matrix Q is defined in (^). Using this, the fermionic path integral 
for Dirac fermions can be written as 



2 

/ [##]e~ 5 ' = / [##]e"^ w = det Q = f[ J [dV^e'^ Q * U) . (35) 

3=1 

For Majorana fields the path integral involves only [d\E^')], either with j — 1 or j — 2. 
For \l> = \l> (1) or ^ = \I> (2) we have 



[dtf] e -2*<2* = iJdetg . (36) 



In order to define the sign in (|36f) one has to consider the Pfaffian of the anti- 
symmetric matrix 

M = CQ. (37) 

This can be defined for a general complex antisymmetric matrix M a p = —Mp a with 
an even number of dimensions (1 < a, [3 < 2N) by a Grassmann integral as 

Pf (M) = J [d<t>]e-^ M ^ = ^e 9lft ... a ^M aift . . . M aN p N . (38) 

Here, of course, [dip] = d<p2N ■ ■ ■ d<f>i, and e is the totally antisymmetric unit tensor. 
One can easily show that 

[Pf(M)] 2 = detM . (39) 

If M is taken from (|37D one also has det M = det Q. 

The relations in (|36D or (^) show that, in order to represent a Majorana fermion, 
in the path integral over the gauge field the square root of the fermion determinant 
(or the Pfaffian of the matrix in (|37D) has to be taken. In this sense a Majorana 
fermion corresponds to a flavour number |. Concerning the sign of the square root, in 
numerical simulations it is easier to take always the absolute value. This presumably 
does not have an influence in the continuum limit because in the continuum the 
(real) eigenvalues of the Dirac matrix come in pairs and the square root is always 
positive (see, for instance, ]T8"||). 

In general, the lattice action describing a given "target" continuum quantum 
field theory is not unique. Besides the Curci-Veneziano action discussed up to now, 



another possibility is based on five- dimensional domain walls |19| , |20| , gT| . In this 
approach one knows the value of the bare fermion mass where the supersymmetric 
continuum limit is best approached and one has advantages from the point of view 
of the speed of symmetry restoration. The price one has to pay is the proliferation of 
(auxiliary) fermion flavours. Another proposal for reaching supersymmetric quantum 
field theories is to try direct dimensional reduction on the lattice |22[ . 



In order to perform Monte Carlo simulations with effective flavour number | cor- 
responding to Majorana fermions in the lattice formulation of Curci and Veneziano, 
one can either use the multi-bosonic technique [p3j |l7f or apply the hybrid classical 



dynamics algorithm [pl |. Exploratory studies have been started recently. (For a 



recent review and status report see |25[ . 
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The first step in numerical simulations is to consider the quenched approximation, 



which neglects the dynamical effects of gluinos |26], |27|. Since quenching breaks su- 
persymmetry explicitly, this mainly serves as a testing ground for mass measurements 
and helps to localize the physically interesting bare parameter range. 

A first large scale numerical simulation of SU(2) SYM with dynamical gluinos has 
been started recently by the DESY-Miinster collaboration [[28[] using the supercom- 



puters at HLRZ, Jiilich and DESY, Zeuthen. The main goals of this collaboration 
are: to find the first order phase transition at zero gluino mass and to determine the 
masses of the lowest bound states formed out of gluons and gluinos in the interesting 
range of the gluino mass. 
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